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ABSTRACT 

We report the detection of very-high-energy (VHE) gamma-ray emission from supernova remnant 
(SNR) G106.3+2.7. Observations performed in 2008 with the VERITAS atmospheric Cherenkov 
gamma-ray telescope resolve extended emission overlapping the elongated radio SNR. The 7.3a (pre- 
trials) detection has a full angular extent of roughly 0.6° by 0.4°. Most notably, the centroid of 
the VHE emission is centered near the peak of the coincident 12 CO (J=l-0) emission, 0.4° away 
from the pulsar PSR J2229+6114, situated at the northern end of the SNR. Evidently the current- 
epoch particles from the pulsar wind nebula are not participating in the gamma-ray production. 
The VHE energy spectrum measured with VERITAS is well characterized by a power-law dN/dE = 
N {E/3 TeV)- T with a differential index of V = 2.29 ± 0.33 stat ± 0.30 sys and a flux of N = (1.15 ± 
0.21 s tat ± 0.35 sas ) x 10~ 13 cm~ 2 s _1 TcV _1 . The integral flux above 1 TeV corresponds to ~ 5 percent 
of the steady Crab Nebula emission above the same energy. We describe the observations and analysis 
of the object and briefly discuss the implications of the detection in a multiwavelength context. 
Subject headings: gamma rays: observations — ISM: individual (G106.3+2.7=VER J2227+608) - 
pulsars: individual (J2229+6114) — supernova remnants 
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1. INTRODUCTION 

The supernova remnant (SNR) G106.3+2.7 was first 
identified as a faint extended radio source in a northern 
Galactic Plane radio survey with the Dominion Radio 
Astrop hysical Observatory (DRAO) by IJoncas fc Higgj 
(1990) . Wor k ing w ith the earlier radio data of 
Kallas & Reich (1980), those authors determined that 
the object was a previously undetected supernova rem- 
nant, with a radio spectral index of 0.45 and a flux 
density at 1 GH z of 6 Jy. Subsequent work by 
IPineault fc Joncasl (|2000D confirmed the object as a su- 
pernova remnant, with an estimated age and distance of 
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1.3 Myr, and 12 kpc. The object can be described as 
cometary in shape, with a compact radio-bright head to 
the north, and a dimmer, extended tail to the southwest. 

Located at the northern edge of the remnant's head 
is the pulsar PSR J2229+6114 and its associated 
boomerang-shape d radio and x-ray-emi tting wind neb- 
ula, G106. 6+2.9 (jHalpern et alj l2001bO . The pulsar, 
which has a period of 51.6 ms, was discovered in a search 
of the error box o f the EG RET source 3EG J2227+6122 
by lHalpern et alJ (|2001al ) and is notable for being one 
of the most energetic pulsars known, with a spin-down 
power E = 2.2 x 10 37 erg s" 1 , d espite having a low radio 
luminosity (jKothes et al.ll2006h . 

Based on polari zation studies and v elocity maps of HI 
and CO emission. iKothes et alJ (|2001h associate the pul- 
sar wind nebula and the supernova remnant with the 
same progenitor event and invoke a scenario of shock- 
wave evolution into a complex molecular environment to 
explain the unusual configuration of the constituent com- 
ponents. These authors derive a distance to the remnant 
of 0.8 kpc, a nd adopt an age of 10 kyr, inferred from the 
pulsar data (jHalpern et al.ll2001af l . 

At GeV energies, the EGRET source 
3EG J2227+6122, with relatively large error con- 
tours, is compatible with the pulsar pos ition, as well as 
with the main bulk of the radio remnant (Har tman et alJ 
1999). The newly-released Large Area Telescope Bright 
Source List of the Fermi Gamma-ray Space Telescope 
(FGST) contains the object OFGL J2229.0+6114, the 
positi on of which is co nsistent with the pulsar posi- 
tion (lAbdo et alJ l2009bf) The AGILE source 1AGL 
J2231+6109 ()Pittori et"atll2009S ) is also located nearby. 

At TeV energies, the most stri ngen t flux limits come 
from the MAGIC Collaboration (|Aliull2007h . who quote 
a point-source upper limit of ~ 10% of the Crab Neb- 
ula flux above 220 GeV at the pulsar position. Less 
constraining point-so urce limits have also been pre- 
sented bv the Whipple (jFegan et al.ll2005h and VERITAS 
(|Konopelko et al.ll2007l ) Collaborations. 

Additionally, the Milagro Collaboration has recently 
reported a new analysis of their data for this region, 
resulting in the detection of gamma rays over a broad 
~ 1° area encompassing both the pulsar positio n and 
the main bulk of the remnant (jAbdo et al.ll2009al >. The 
detection corresponds to the loca tion of their previ ously 
published "source candidate C4" (jAbdo et al.ll2007h . but 
lacks the angular resolution to provide a definitive asso- 
ciation with a particular region within the SNR/pulsar 
complex. The statistical significance of the detection at 
the pulsar location is 6.6cr and the median energy of the 
detected gamma rays is 35 TeV. The quoted differential 
flux is (70.9±10.8 stat ±35.5 sys ) x 10~ 17 TeV" 1 cm- 2 s _1 , 
though this depends on assumptions about the source 
spectral shape. 

2. VERITAS INSTRUMENT & OBSERVATIONS 

The VERITAS observatory comprises an array of four 
12 m imaging atmospheric Cherenkov telescopes, located 
at Mt. Hopkins in southern Arizona (1268 m a.s.l., 
N 31°40'30", W 110°57'07"). The instrument, which 
features a 3.5° field-of-view, has a 5cr point source sen- 
sitivity of 1% of the steady Crab Nebula flux in under 
50 hours and an angular resolution of ~ 0.1°, making it 
well-suited for the study of extended gamma-ray sources. 



For a detai led description of the tel escopes and t h eir op - 
eration see iWeekes et all (|200l or Holder et~aLl ([2006). 
A recent review of the atmospher i c Cher enkov technique 
can be found in lAharonian et al.1 (|2008ch . 

VERITAS observed G106.3+2.7 during three epochs. 
The first was during November/December of 2006, with 
2 telescopes , resulting in a point-so urce flux upper limit 
reported in iKonopelko et all (|2007l ). Observations with 
4 telescopes continued in September-November of 2007, 
as a part of a larger survey for PWN emission from high- 
spin- powered pulsars in the northern sky (jAliu et alJ 
2008|). After a hmt of extended gamma-ray emission was 
seen in these data, follow-up observations were triggered 
in the 2008 observing season. In each epoch, the tele- 
scope tracking direction was centered on the position of 
the pulsar. 

All observations were performed in the so-called "wob- 
ble mode" , where the telescopes point in a direction offset 
by some angle from the putative source position. This is 
done to provide a simultaneous estim ate of backgrounds 
with every run (|Fomin et alJ [T994). The offset direc- 
tion is sequentially varied from run to run between the 
four cardinal directions. In the earlier observing epochs, 
small offset angles, suitable for point-source searches, 
were used. Based on the traces of extended emission seen 
in the previous data, a 0.7 degree offset, more appropriate 
for larger sources, was adopted in the 2008 epoch. Only 
data from this epoch, which also includes the deepest 
exposure, are included in this Letter. After quality se- 
lection cuts to remove periods with poor weather or mal- 
functioning hardware, 111 observing runs remain, with 
a mean zenith angle of 34 degrees and a total deadtimc- 
corrected exposure of 33.4 hours. 

3. DATA ANALYSIS AND RESULTS 

The data were analyzed using t he VERITAS stan - 
dard analysis package, d escribed in iCogan et alJ (|2007l ) 
and lAcciari et al.l (|2008l ). After gain corrections and 
cleaning, the air shower images in the individual tele- 
scop es are param eterized using the Hillas moment anal- 
ysis (|Hillasl 119851 ). Images which pass the pre-selection 
criteria (number of focal-plane phototubes in image 5; 
distance from image centroid to camera center < 1.43 
degrees; total image size > 150 photoelectrons) are then 
retained for full event reconstruction using the mean- 
scaled width and len g th pa rameters (see e.g., [Konopclko 
(1995); iDaniel et all (|2007l )). Events where only a sin- 
gle telescope image survives cuts are rejected, as are 
events where the only surviving images are from tele- 
scopes Tl and T4, which, with a separation of 35 m, are 
too close together for reliable event reconstruction. The 
gamma/hadron cuts were based on an optimization us- 
ing a large sample of Crab Nebula data, scaled to a flux 
level of 3%. The energy threshold of the analysis is 630 
GeV. 

Background estimation for the spectral reconstruc- 
tion was performed u sing the reflected-region model 
(|Aharonian et al.l 120 01). For sky map construction (see 
Figure [l}, the ring- background model was employed 
(for a discussion, see iBerge et alJ (|2007l ) and references 
therein). In the search for extended emission, a squared 
angular integration radius of 0.08 deg 2 , well-suited for 
the investigation of larger sources, was used. Addition- 
ally, a smaller radius was used in a point source search. 
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Fig. 1. — Sky map of TeV gamma-ray emission from G106.3+2.7, as measured by VERITAS. The color scale indicates the number of 
excess gamma-ray events from the region, using a squared integration radius of 0.08 deg 2 . The centroid of the TeV emission is indicated wit h 
a thin black cross. Overlaid are 1420 MHz radio contours from the DRAO Synthesis Telescope (thin black lines - Landcckcr ct al. (2000)) 
and 12 CO emission (J=l-0) from the high-resolution FCRAO Survey (magenta lines - -Hey er e t al., (1993)). The CO velocity selection is 
discussed in the text. The open yellow cross shows the location of pulsar PSR J2229+6114 and the yellow star at the left is the AGILE 
bright source 1AGL J2231+6109. The yellow circle indicates the 95% error contour for the Fermi source 0FGL J2229. 0+6114. The circle 
labeled PSF represents the VERITAS gamma-ray point-spread function for this analysis (68% containment). 



The statistical significance of any exc ess event counts 
were evaluated using Equation (17) from lLi fc Mai (fl983h 
and the results were confirmed using an independent in- 
ternal analysis chain. There are no bright stars in the 
field of view to complicate observations or analysis. 

Using the extended emission cuts, a pre-trials maxi- 
mum significance of 7.3tr was observed. The location of 
this peak significance is ~ 0.4 degrees to the south of the 
position of the pulsar PSR J2229+6114. Because there is 
no a priori test point at this location, we apply a conser- 
vative trials penalty determined by tiling the entire area 
defined by the remnant's outer radio contour with 0.04° 
square search bins. An additional factor associated with 
the small set of search cuts is also incorporated. With 
these included, the post-trials significance of the peak 
emission point is 6.0a. The post-trials significance of the 
excess seen at the nominal pulsar location is 3.6a. 



3.1. Morphology 



Inspection of Figure [T] shows that the angular extent 
of the TeV emission region exceeds the size of the 0.11° 
VERITAS point-spread-function (PSF; the radius that 
contains 68% of the events coming from a point source). 
While the shape of the extension in this map is ap- 
parently more complex than a simple Gaussian in na- 
ture, a statistically acceptable fit can nevertheless be ob- 
tained by fitting a 2-dimensional Gaussian to the uncor- 
rected acceptance-corrected map of excess event counts, 
binned in 0.05° bins. The result of this analysis shows 
that the TeV emission, after accounting for the PSF 
of the instrument, can be characterized by a la angu- 
lar extent of 0.27° ± 0.05° along the major axis, and 
0.18° ± 0.03° along the minor axis; the orientation an- 
gle is 22° east of north. The centroid of the fit lies at 
22 h 27 m 59 s , +60° 52' 37" (J2000) and hence we assign 
the identifier VER J2227+608. The statistical uncer- 
tainty in the centroid position is 0.07° in RA and 0.04° 
in declination, with an additional combined systematic 
uncertainty of 0.07°. 
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3.2. Spectrum 

The spectrum of the source was determined by inte- 
grating over a circular region of radius 0.32°, positioned 
to encompass the majority of the observed emission. The 
resulting spectrum, which is shown in Figure can be 
fitted (% 2 =0.76;ndf=3) with a simple power-law of the 
form dN/dE = N {E/3 TeV)~ T ', with a differential flux 
constant of N = (1.15 ± 0.27 stat ± 0.35 sys ) x 10~ 13 
crxi- 2 s- 1 TeV" 1 and an index of T = 2.29 ± 0.33 stat ± 
0.30^. The total flux integrated over the remnant for 
energies above 1 TeV is (l.ll±0.25 stat ±0.28 s2/s ) x 1(T 12 
cm -2 s -1 , which corresponds to ~ 5% of the steady Crab 
Nebula flux . Also shown in this figure is the Milagro flux 
point from lAbdo et al.1 (|2009al ). The energy resolution 
of this analysis is better than ~ 20% for energies above 
1 TeV. 

4. DISCUSSION 

Figure Q] shows the TeV gamma-ray image for 
G106.3+2.7. The color scale indicates the number of ex- 
cess gamma-ray candidate events, as measured by VER- 
ITAS. Radio continuum emission at 1420 MHz mea- 
sured by the Dominion Radio As trophysical Obse r vator y 
(DRAO) Synthesis Telescope (|Landecker et al.l 120001) . 
and supplied through t he Canadian Galactic Plane Sur- 
vey (| Taylor et al. 2003), is indicated with thin black con- 
tours. A high-resolution map of 115 GHz line emission 
associated with the molecular 12 CO (J=l-0) transition, 
from the Five College Radio Astronomy Observatory 
(FCR AO) is shown in magenta contours dHever et al.l 
119981) . Following the analysis of iKothes et all (|2001f ), 
the CO velocities in the map are restricted to between 
-4 km s _1 and -6 km s _1 . This slice encompasses the ma- 
jority of the cloud emission overlapping the radio rem- 
nant. Note that the CO contours have been smoothed 
with a ~ 3' Gaussian kernel. The open yellow cross 
marks the location of PSR J2229+6114; its wind neb- 
ula, the "Boomerang" (~ 0.08° in diameter), can be 
seen in the radio contours around it. The yellow star at 
the left is the AGILE bright source 1AGL J2231+6109 
and the 95% FGST confidence contour for the source 
0FGL J2229.0+6114 is shown as the yellow circle cen- 
tered on the pulsar location. 

As can be seen from Figure [1] the TeV emission ap- 
pears correlated with the position of the radio contours 
of the main body of the remnant, reaching maximum 
flux ~ 0.4° away from the pulsar. Indeed, the region of 
greatest intensity is apparently coincident with the loca- 
tion of the molecular cloud emission. This coincidence 
is re miniscent of b e havio r seen in similar o b jects ( see, 
e.g., lAcciari et all (|2009f k lAharonian et al.l (|2008al lbl)) 
and is supportive of theories of had ronically-induced 
(that is, p + p — » 77° — > 77) emi ssion (jAharonian et al.l 
119941 lAharonian fc Atovan|[l99l ). The detection of OH 
maser emission, a signpost of shock wave/cloud interac- 
tions dElitzuHU976tlWardle fc Yusef-Zadehll2002D . would 
strengthen this claim, but no such observations have been 
reported. 

A leptonic scenario, wherein gamma rays are pro- 
duced by the inverse-Compton scattering of a local pho- 
ton pool (for instance, synchrotron-produced photons, 
cosmic microwave background (CMB) photons, or in- 
frared photons) by high-energy electrons, produced ei- 
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Fig. 2. — Differential energy spectrum of G106.3+2.7 as measured 
by VERITAS. The error bars represent lcr statistical errors only. 
The solid line shows the results of a pure power-law fit and the 
dotted line shows the flux of the Crab Nebula for comparison. Als o 
shown (triangle) is the Milagro flux point from Abdo et al. (2009a), 
plotted with summed systematic and statistical errors. The lower 
panel shows the residuals of the VERITAS data to the power-law 
fit. The details of the analysis and the fit are discussed in the text. 

ther by the SNR, or the P WN, is another possibility. 
iBednarek fc BartosiU (|2005f ) have applied such a model 
to the Boomerang PWN and made predictions for the 
expected gamma-ray flux. In this model, leptons can ob- 
tain high energies through resonant scattering off Alfven 
waves generated by heavy nuclei, which, in turn, have 
been accelerated directly by the pulsar or PWN. Assum- 
ing only photon targets of synchrotron and CMB ori- 
gin, they obtain an integral flux prediction which cor- 
responds to (in the absence of any cutoff) 1.6 x 10~ 12 
photons cm~ 2 s _1 above 1 TeV. This prediction is con- 
sistent within statistical and systematic errors with the 
present measurement. The lack of a clear spatial co- 
incidence between the PWN and t he TeV emission is 
not unique {e.g., lAharonian et alj (2006a. b)) and can 
be explained in an older (r > 10 4 yrs) remnant such 
as this because the emission is dominated by the scat- 
tering of CMB photons by lower-energy leptons. These 
leptons may be displaced from the present pulsar posi- 
tion due either to pulsar m otion or to interactions be- 
tween the PWN and SNR (jBednarek fc Bartosikl 120031 
IBartko fc Bednarekj[2008l) . A possible gamma-ray com- 
ponent due to hadrons trapped i n nearby molecular 
clouds could also contribute {e.g., IBartko fc Bednarekl 
(2008)). Finally, we note that our measured spectrum, 
as shown in Figure if extrapolated to 35 TeV, is con- 
sistent within errors with t h e flux reported by the Mi- 
lagro group in lAbdo et alj (|2009aD . Given the uncer- 
tainty in the Milagro source position, it isn't possible 
to unambiguously associate the two sources. However, if 
there is, in fact, a single mechanism at work in producing 
the emission seen by both instruments, then this could 
weaken the case for a leptonic scenario, given that the 
combined spectrum lacks the curvature expected from 
a ty pical inverse Compton model at these energies (see, 
e.q. jde Jager fc Diannati-Ataj (|2008fl ). 

5. CONCLUSION 

The VERITAS observatory has made a high-confidence 
detection of TeV gamma rays from the supernova rem- 
nant G106.3+2.7. This gamma-ray emission is signifi- 
cantly extended, and seems well-correlated with the ra- 
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dio extent of the remnant, in particular, with the location 
of known nearby molecular clouds. The emission seen at 
the location of the associated pulsar PSR J2229+6114 is 
non-zero, but only of marginal statistical significance. 

The total flux from the remnant above 1 TeV is about 
~ 5% of the Crab Nebula, and can be described by a 
power-law in energy with differential index 2.29. With 
the information available, it is not yet possible to make a 
firm determination of whether the TeV gamma-ray emis- 
sion is of leptonic or hadronic origin, though a solid as- 
sociation with the Milagro source could favor the latter. 
With additional exposure and a refined observation strat- 
egy (i.e., wobbling around the location of the emission 
peak rather than the pulsar) a more in-depth study of 
the object will be achievable, including an investigation 
of possible energy dependence in the source morphology. 
Such an effort will benefit from the enhanced sensitivity 
and angular resolution of the reconfigured VERITAS ar- 



ray, in which the Tl telescope (see Section[3|) is relocated 
to obtain a more optimal baseline. This upgrade will be 
completed in autumn, 2009. 
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